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Abstract 
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on the latter, as introduced by Crnkovic-Witten and Zuckerman, from the multi- 
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symplectic structure, according to the standard rules, is precisely the covariant 
bracket due to Peierls and DeWitt. 
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1 Introduction 



One of the most annoying flaws of the usual canonical formalism in field theory is 
its lack of manifest covariance, that is, its lack of explicit Lorentz invariance (in the 
context of special relativity) and more generally its lack of explicit invariance under 
space-time coordinate transformations (in the context of general relativity). Of course, 
this defect is built into the theory from the very beginning, since the usual canonical 
formalism represents the dynamical variables of classical field theory by functions on 
some spacelike hypersurface (Cauchy data) and provides differential equations for their 
time evolution off this hypersurface: thus it presupposes a splitting of space-time into 
space and time, in the form of a foliation of space-time into Cauchy surfaces. As a result, 
canonical quantization leads to models of quantum field theory whose covariance is far 
from obvious and in fact constitutes a formidable problem: as a well known example, we 
may quote the efforts necessary to check Lorentz invariance in (perturbative) quantum 
electrodynamics in the Coulomb gauge. 

These and similar observations have over many decades nourished attempts to de- 
velop a fully covariant formulation of the canonical formalism in classical field theory, 
which would hopefully serve as a starting point for alternative methods of quantization. 
Among the many ideas that have been proposed in this direction, two have come to 
occupy a special role. One of these is the "covariant functional formalism" , based on 
the concept of "covariant phase space" which is defined as the (infinite-dimensional) 
space of solutions of the equations of motion. This approach was strongly advocated 
in the 1980's by Crnkovic, Witten and Zuckerman [1-3] (see also [4]) who showed how 
to construct a symplectic structure on the covariant phase space of many important 
models of field theory (including gauge theories and general relativity), but the idea 
as such has a much longer history. The other has become known as the "multisymplectic 
formalism" , based on the concept of "multiphase space" which is a (finite-dimensional) 
space that can be defined locally by associating to each coordinate q % not just one 
conjugate momentum pi but n conjugate momenta pf (/i — 1, . . . , n), where n is the 
dimension of the underlying space-time manifold. In coordinate form, this construction 
goes back to the classical work of De Donder and Weyl in the 1930's [5,6], whereas 
a global formulation was initiated in the 1970's by a group of mathematical physi- 
cists, mainly in Poland [7-9] but also elsewhere [10-12], and definitely established in 
the 1990's [13,14]; a detailed exposition, with lots of examples, can be found in the 
GIMmsy paper [15]. 

The two formalisms, although both fully covariant and directed towards the same 
ultimate goal, are of different nature; each of them has its own merits and drawbacks. 

• The multisymplectic formalism is manifestly consistent with the basic principles of 
field theory, preserving full covariance, and it is mathematically rigorous because 
it uses well established methods from calculus on finite-dimensional manifolds. 
On the other hand, it does not seem to permit any obvious definition of the 
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Poisson bracket between observables. Even the question of what mathematical 
objects should represent physical observables is not totally clear and has in fact 
been the subject of much debate in the literature. Moreover, the introduction 
of n conjugate momenta for each coordinate obscures the usual duality between 
canonically conjugate variables (such as momenta and positions), which plays a 
fundamental role in all known methods of quantization. A definite solution to 
these problems has yet to be found. 

• The covariant functional formalism fits neatly into the philosophy underlying the 
symplectic formalism in general; in particular, it admits a natural definition of 
the Poisson bracket (due to Peierls [16] and further elaborated by DeWitt [17-19]) 
that preserves the duality between canonically conjugate variables. Its main draw- 
back is the lack of mathematical rigor, since it is often restricted to the formal 
extrapolation of techniques from ordinary calculus on manifolds to the infinite- 
dimensional setting: transforming such formal results into mathematical theorems 
is a separate problem, often highly complex and difficult. 

Of course, the two approaches are closely related, and this relation has been an 
important source of motivation in the early days of the theory [8]. Unfortunately, 
however, the tradition of developing them in parallel seems to have partly fallen into 
oblivion in recent years, during which important progress was made in other directions. 

The present paper, based on the PhD thesis of the second author [21], is intended to 
revitalize this tradition by systematizing and further developing the link between the 
two approaches, thus contributing to integrate them into one common picture. It is or- 
ganized into two main sections. In Sect. 2, we briefly review some salient features of the 
multisymplectic approach to geometric field theory, focussing on the concepts needed 
to make contact with the covariant functional approach. In particular, this requires a 
digression on jet bundles of first and second order as well as on the definition of both 
extended and ordinary multiphase space as the twisted affine dual of the first order 
jet bundle and the twisted linear dual of the linear first order jet bundle, respectively: 
this will enable us to give a global definition of the space of solutions of the equa- 
tions of motion, both in the Lagrangian and Hamiltonian formulation, in terms of a 
globally defined Euler - Lagrange operator £ and a globally defined De Donder - Weyl 
operator D, respectively. To describe the formal tangent space to this space of solu- 
tions at a given point, we also write down the linearization of each of these operators 
around a given solution. In Sect. 3, we apply these constructions to derive a general 
expression for the symplectic form Q on covariant phase space, a la Crnkovic-Witten- 
Zuckerman, in terms of the multisymplectic form uj on extended multiphase space. 
Then we prove, as the main result of this paper, that the Poisson bracket associated 
with the form Q, according to the standard rules of symplectic geometry, suitably ex- 
tended to this infinite-dimensional setting, is precisely the Peierls- DeWitt bracket of 
classical field theory [16-19]. Finally, in Sect. 4, we comment on the relation of our 
results to previous work and on perspectives for future research in this area. 
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2 Multisymplectic Approach 



2.1 Overview 



The multisymplectic approach to geometric field theory, whose origins can be traced 
back to the early work of Hermann Weyl on the calculus of variations [6], is based on 
the idea of modifying the transition from the Lagrangian to the Hamiltonian framework 
by treating spatial derivatives and time derivatives of fields on an equal footing. Thus 
one associates to each field component ip 1 not just its standard canonically conjugate 
momentum n i but rather n conjugate momenta irf , where n is the dimension of space- 
time. In a first order Lagrangian formalism, where one starts out from a Lagrangian L 
depending on the field and its first partial derivatives, these are obtained by a covariant 
analogue of the Legendre transformation 

This allows to rewrite the standard Euler- Lagrange equations of field theory 

o 9L dL 

as a covariant first order system, the covariant Hamiltonian equations or De Donder- 
Weyl equations 

dH j dH 

tot= 9 »<P ' dji = ( 3 ) 

where 

H = nfd^ - L (4) 

is the covariant Hamiltonian density or De Donder- Weyl Hamiltonian. 

Multiphase space (ordinary as well as extended) is the geometric environment built 
by appropriately patching together local coordinate systems of the form {q l ,Pi) - in- 
stead of the canonically conjugate variables (q l ,Pi) of mechanics - together with space- 
time coordinates and, in the extended version, a further energy type variable that 
we shall denote by p (without any index). The global construction of these multi- 
phase spaces, however, has only gradually come to light; it is based on the following 
mathematical concepts. 



• The collection of all fields in a given theory, defined over a fixed (n-dimensional 
orientable) space-time manifold M, is represented by the sections ip of a given 
fiber bundle E over M, with bundle projection ir : E — > M and typical fiber Q. 
This bundle will be referred to as the configuration bundle of the theory since Q 
corresponds to the configuration space of possible field values. 
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• The collection of all fields together with their partial derivatives up to a certain 
order, say order r, is represented by the r-jets j r ip = (ip, dp, . . . , d r ip) of sections 
of E, which are themselves sections of the r th order jet bundle J r E of E, regarded 
as a fiber bundle over M. In this paper, we shall only need first order jet bundles, 
with one notable exception: the global formulation of the Euler - Lagrange equa- 
tions requires introducing the second order jet bundle. 

• Dualization - the concept needed to pass from the Lagrangian to the Hamiltonian 
framework via the Legendre transformation - comes in two variants, based on the 
fundamental observation that the first order jet bundle J l E of E is an affine 
bundle over E whose difference vector bundle J 1 E will be referred to as the 
linear jet bundle. Ordinary multiphase space is obtained as the twisted linear 
dual J 1@ E of J 1 E while extended multiphase space is obtained as the twisted 
affine dual J l ®E of J l E, where the prefix "twisted" refers to the necessity of 
taking an additional tensor product with the bundle of n-forms on M. 1 

• The Lagrangian £ is a function on J l E with values in the bundle of n-forms 
on M so that it may be integrated to provide an action functional which enters 
the variational principle. The De Donder-Weyl Hamiltonian Oi is a section of 
J l ®E, considered as an affine line bundle over J 1 ®E. 

Note that the formalism is set up so as to require no additional structure on the con- 
figuration bundle or on any other bundle constructed from it: all are merely fiber 
bundles over the space-time manifold M. Of course, additional structures do arise 
when one is dealing with special classes of fields (matter fields and the metric tensor 
in general relativity are sections of vector bundles, connections are sections of affine 
bundles, nonlinear fields such as those arising in the sigma model are sections of trivial 
fiber bundles with a fixed Riemannian metric on the fibers, etc.), but such additional 
structures depend on the kind of theory considered and thus are not universal. Finally, 
the restriction imposed on the order of the jet bundles considered reflects the fact 
that almost all known examples of field theories are governed by second order partial 
differential equations which can be derived from a Lagrangian that depends only on the 
fields and their partial derivatives of first order, which is why it is reasonable to develop 
the general theory on the basis of a first order formalism, as is done in mechanics [22,23]. 

2.2 The First Order Jet Bundle 

The field theoretical analogue of the tangent bundle of mechanics is the first order jet 
bundle J l E associated with the configuration bundle E over M. Given a point e in E 

1 We use an asterisk * to denote linear duals of vector spaces or bundles and a star ★ to denote affine 
duals of affine spaces or bundles. These symbols are appropriately encircled to characterize twisted 
duals, as opposed to the ordinary duals defined in terms of linear or affine maps with values in R. 
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with base point x = n(e) in M, the fiber J\E of J l E at e consists of all linear maps 
from the tangent space T X M of the base space M at x to the tangent space T e E of the 
total space E at e whose composition with the tangent map T e ir : T e E — * T^M to the 
projection 7r : E — > M gives the identity on T^M: 

= { 7 eL(T x M,T e E) : T>o 7 = id T ^ M } . (5) 

Thus the elements of JlE are precisely the candidates for the tangent maps at x to 
(local) sections (p of the bundle E satisfying (p(x) = e. Obviously, J\E is an afline 
subspace of the vector space L(T X M, T e E) of all linear maps from T X M to the tangent 
space T e E, the corresponding difference vector space being the vector space of all linear 
maps from T X M to the vertical subspace V e E: 

JlE = { 7 G L(T X M, T e E) : 7> o 7 = 0} = L(T X M, V e E) = T* X M ®V e E . (6) 

The jet bundle J l E thus defined admits two different projections, namely the target 
projection r E : J 1 E — > E and the source projection a E : J 1 E — > M which is simply 
its composition with the original bundle projection, that is, a E = n o r E . The same 

— » . 

goes for J E, which we shall call the linearized first order jet bundle or simply linear 
jet bundle associated with the configuration bundle E over M. 

The structure of J l E and of J l E as fiber bundles over M with respect to the source 
projection (in general without any additional structure), as well as that of J 1 E as an 
afline bundle and of J 1 E as a vector bundle over E with respect to the target projection, 
can most easily be seen in terms of local coordinates. Namely, local coordinates x^ for M 
and q % for Q, together with a local trivialization of E, induce local coordinates (x M , q l ) 
for E as well as local coordinates {x 11 , q l , q 1 ^) for J l E cL(tt*(TM),TE) and (x^,q\ q^) 
for J l E c L(n*(TM),TE). Moreover, local coordinate transformations x^ — > x' v for M 
and q l — > g /J for Q, together with a change of local trivialization of E, correspond to 
a local coordinate transformation (x M , g l ) — > (x /iy , g /J ) for E where 

x» = x' u {x^) , q' j = q' j {x^,q l ) . (7) 

The induced local coordinate transformations (x^, q l , q l ) — > {x' u ,q'i ,q'J) for J 1 ^ and 
(x M , (f , qp — ■> (x /iy , g'- 7 , g^) for J 1 ^ are then easily seen to be given by 

fj dx» dq'i t dx» dq'i 
qv dx>» 8$ % dx'» dx» ' K ' 



and 



dx^dq^ 
'" dx'» ()(/' • 
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This makes it clear that J l E is an affine bundle over E with difference vector bundle 



J l E = T*M ® VE , (10) 

in accordance with eq. 

That the (first order) jet bundle of a fiber bundle is the adequate arena to incorpo- 
rate (first order) derivatives of fields becomes apparent by noting that a global section 
ip of E over M naturally induces a global section j l (p of J l E over M given by 

j l <p(x) = 7> 6 Jl [x) E for xeM . 

In the mathematical literature, j 1 ^ is called the (first) prolongation of ip, but it would 
be more intuitive to simply call it the derivative of ip since in the local coordinates used 
above, 

Mx) = (a* ¥>«(*), fl^x)) , 
where <9 M = d/dx^; this is symbolically summarized by writing j l p = (<p,dip). 

Similarly, it can be shown that the linear jet bundle of a fiber bundle is the adequate 
arena to incorporate covariant derivatives of sections, with respect to an arbitrarily 
chosen connection. 

Finally, let us discuss briefly the lifting, from E to J X E, of (local) bundle automor- 
phisms and, passing to generators of one-parameter groups, of projectable vector fields. 
Let $ : E —>■ E be an automorphism of the fiber bundle E over M and <p : M — > M 
the induced diffeomorphism of M such that the diagram 

$ 

E — > E 

7T 7T 

M — > M 



commutes. This can be lifted to an automorphism of the jet bundle J X E, as an affine 
bundle over E, by defining J 1( 3> : J X E — > J l E as follows: given a point e in E with 
base point x = 7r(e) in M and a 1-jet 76 J^E, define the 1-jet J x $(7) e J\,^E by 

J^irt) = T e $ o 7 o {T x 4>Y l . (11) 

2 Given any vector bundle V over M, such as TM, T*M or any of their exterior powers, one can 
consider it as as vector bundle over E by forming its pull-back tt*V. In order not to overload the 
notation, we shall here and in what follows suppress the symbol ir*. 
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Obviously, this formula defines a linear map from L(T x M,T e E) to L(T^ x ^M,T^ e -jE) 
that restricts to an affine map from J\E to J\t e \E : 



T* (e )7T O J X $( 7 ) 



T*( e )iroT e $ 070 (T x .0) 

T e (7TO$)o 7 o(T»- 1 
T e (0O7T) 7 (T» _1 

T x (p o T e 7i o ^ o {T^y 1 



id 



0(a) 



In particular, the diagram 



J 1 ® 



J X E 



J l E 



T 



E 



E 



commutes, justifying to call J 1( 3> the first prolongation of $. This construction can 
be generalized to any fiber bundle map $ : E — > F over a (local) diffeomorphism 
: M — > N , giving an affine bundle map J x $ : J 1 ^ — > J 1 / 71 over § : E F. 

Passing to the description of the infinitesimal situation, let us consider a projectable 
vector field V on E 1 , whose flow is a one-parameter group of (local) automorphisms of E 
that can be lifted to a one-parameter group of local automorphisms of J l E, generated 
by a projectable vector field J l V on J X E: this is then defined to be the prolongation 
of V. Thus 



V 



9$, 



jV 



d{J 1 ^? > 



A=0 



<9A 



A=0 



In local coordinates as before, we can write 



d ■ d 

V = V + V* r 

dx^ dq 1 



12) 



where = V tl (x u ) and V % = V l {x v , qi), and since the lifting of bundle automorphisms 
is described by the transformation law (JEJ), differentiation with respect to A gives 



ox^ oq l 



dV 



dq 



k 



dV K 
dx^ 



dV<\ d 



dx^ J dq 



(13) 



2.3 Duality 



The next problem to be addressed is how to define an adequate notion of dual for J X E. 
The necessary background information from the theory of affine spaces and of affine 
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bundles (including the definition of the affine dual of an affine space and of the transpose 
of an affine map between affine spaces) is summarized in the Appendix. Briefly, the 
rules state that if A is an affine space of dimension k over R, its dual A* is the space 
A(A,M) of affine maps from A to R, which is a vector space of dimension k + 1. 
Thus the affine dual J U E of J l E and the linear dual J l *E of J l E are obtained by 
taking their fiber over any point e in E to be the vector space 

Jl*E = { z e : JlE — ► R affine } (14) 

and 

J\*E = { z e : J\E — ► R linear} (15) 

respectively. However, as mentioned before, the multiphase spaces of field theory are 
defined with an additional twist, replacing the real line by the one-dimensional space 
of volume forms on the base manifold M at the appropriate point. In other words, the 
twisted affine dual 

ji® E = j^e® f\ n T*M (16) 
of J 1 E and the twisted linear dual 

fi® E = p*E®f\ n T*M (17) 

of J l E are defined by taking their fiber over any point e in E with base point x = 7r(e) 
in M to be the vector space 

Jl®E = {z e : j\E — > A" KM affine} (18) 

and 

J l ®E = {z e :JlE^I\ n T* x M linear} (19) 

respectively. 2 As in the case of the jet bundle and the linear jet bundle, all these 
duals admit two different projections, namely a target projection onto E and a source 
projection onto M which is simply its composition with the original projection n. 

Using local coordinates as before, it is easily shown that all these duals are fiber 
bundles over M with respect to the source projection (in general without any additional 
structure) and are vector bundles over E with respect to the target projection. Namely, 
introducing local coordinates (x M , q l ) for E together with the induced local coordinates 
(x M , q 1 , q l ) for J l E and (a^,(f,g^) for J 1 E as before, we obtain local coordinates 
(x^,q l ,Pi,p) both for J X *E and for J l ®E as well as local coordinates (x M ,g*,pf) both 
for J 1 *E and for J 1 ®E, respectively. These are defined by requiring the dual pairing 
between a point in J l *E or J l ®E with coordinates (x M ,g*,pf,p) and a point in J 1 E 
with coordinates (x M , q\q 1 ^) to be given by 

PH + P ( 2 °) 
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in the ordinary (untwisted) case and by 



(p?<£ + p) d n x (21) 

in the twisted case, whereas the dual pairing between a point in J 1 *E or in J X ®E with 
coordinates (x^,q l ,p^) and a point in J 1 E with coordinates (f, qV) is given by 

pX (22) 

in the ordinary (untwisted) case and by 

Ptq; d n x (23) 

in the twisted case. Moreover, a local coordinate transformation (x^, q % ) — *> (x fu , g /jf ) 
for E as in eq. (JJJ) induces local coordinate transformations for J X E and for J l E as in 
eqs (jHJ) and (JOJ) which in turn induce local coordinate transformations (x^ , q l , pf , p) — > 
(x /iy , q'^ iP'j i p') both for J 1 *^ and for J X ®E as well as local coordinate transformations 
(x (J, ,q 1, ,Pi) — > {x' v ,q'i .,p'j) both for J x *i? and for J l ®E: these are given by 

„, _ cte^ d(f^ jj_ , _ _ chp_ d(f_ „ 

P i ~ dx» dq'o Pl ' P ~ P dx» dq'i Pi [ ] 

in the ordinary (untwisted) case and 

p ? = det (£) d ~h w>** ' pl = det (£) { p -^w p *) (25) 

in the twisted case. 

Finally, it is worth noting that the affine duals J l *E and J X ®E of J l E contain line 
subbundles Jl*E and J^®E whose fiber over any point e in E with base point x = n(e) 
in M consists of the constant (rather than affine) maps from J\E to R and to A" T*M, 
respectively, and the corresponding quotient vector bundles over E can be naturally 
identified with the respective linear duals J X *E and J X ®E of J 1 E, i.e., we have 

J x *E/Jl*E ^ P*E (26) 

and 

J l ®E/J l c ®E S P®E (27) 

respectively. This shows that, in both cases, the corresponding projection onto the 
quotient amounts to "forgetting the additional energy variable" since it takes a point 
with coordinates (rr M , q\p^ ,p) to the point with coordinates (x^,q\pf); it will be 
denoted by r] and is easily seen to turn J X *E and J l ®E into affine line bundles over 
J X *E and over J 1 ®E 1 respectively. 
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2.4 The Second Order Jet Bundle 



For an appropriate global formulation of the standard Euler - Lagrange equations of field 
theory, which are second order partial differential equations, it is useful to introduce 
the second order jet bundle J 2 E associated with the configuration bundle E over M. 
It can be defined either directly, as is usually done, or by invoking an iterative procedure, 
which is the method we shall follow here. Starting out from the first order jet bundle 
J X E of E, regarded as a fiber bundle over M, we consider its first order jet bundle 
J X J X E and define, in a first step, the semiholonomic second order jet bundle J 2 E of E 
to be the subbundle of J X J X E given by 

J 2 E = {neJ l J l E : tji e (k) = J x t e (k)} (28) 

where Tji E : J X J X E — ► J X E is the target projection of J X J X E while J x r E : J X J X E — > 
J X E is the prolongation of the target projection r E : J X E — > E of J X E, considered as 
a map of fiber bundles over M. As will become clear below, J 2 E is an affine bundle 
over E, with difference vector bundle (T*M © (T*M <g> T*M)) <g> VE. Therefore, using 
the construction of the affine quotient of an affine space (bundle) by a vector subspace 
(subbundle) of its difference vector space (bundle), as explained in the Appendix, we 
may complete the construction by observing that since T*M © T*M contains A T*M 
as a vector subbundle (and hence so does T* M © (T* M ©T* M)) , it is possible to define 
the second order jet bundle J 2 E of E as the quotient 

J 2 E = J 2 E I /\ 2 T*M ®VE . (29) 

Once again, J 2 E is an affine bundle over E, with difference vector bundle 

pE = (T*M © V 2 T*M) © VE . (30) 



These assertions can be proved by introducing local coordinates (x^, q l ) for E to- 
gether with the induced local coordinates (x^ 1 , q\ qV) for J X E as before to first define 
induced local coordinates (x^, q' l ,qi, r p , q pp ) for J X J X E. Simple calculations then show 
that the points of J 2 E are characterized by the condition q p = r l p and the points 
of J 2 E by the additional condition q pp = q pp . Moreover, a local coordinate trans- 
formation (x^, q % ) — > (x' u , q'i) for E as in eq. (J7J) induces a local coordinate transfor- 
mation for J X E as in eq. (JHJ) which in turn induces a local coordinate transformation 
(x^, q\ q l , r l , q pp ) — > (x' u , g /jf , q'J , r'J, q'J a ) for J X J X E, given by eq. (JHJ) together with 

„ dx» d<fj , dx* dq'i 

dx' v dq i " dx* ^ ' 1 J 

• dg^ { dx?_ dql 

dx>° dqi q » p dx>° dxf ' 1 ' 
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In particular, eqs (JBJ) and (}3*T]) show that q^ = implies q'J = r'J, as required by 
the global, coordinate independent nature of the definition of J 2 E as a subbundle of 
J X J X E, while eq. (JH2D can be further evaluated by differentiating eq. (JBJ with respect 
to q l and x p , which leads to 

,• dx p dx p dq' j i 
%a = ~dx^ ~dx^ ~dq 1 ?w 

/ dx^ dx p d 2 q' J dx K dx p dx p d 2 x' x dq'^\ i 

+ \ dx' v dx" 7 dx p dq i dx' v dx" 7 dx' x dx p dx K dq i J ?Al 



+ 



<9x p <9 2 g /J dx K dx p dx p d 2 x' x dq 



dx w dx" 7 dx p dx p dx w dx" 7 dx' x dx p dx K dx p 



This is also the induced local coordinate transformation for J 2 E, whereas that for J 2 E 
is obtained by symmetrization: 

<9x M dx p dq' J i 
%a = dx^ dx^ W %p 



1 / dx p dx p dx p dx p \ d 2 q' J dx K dx p dx p d 2 x' x dq 



+ V2 \dx' u dx" 7 + dx" 7 dx' v ) dx p dq 1 dx w dx" 7 dx' x dx p dx K dq 1 J ^ 
/ dx p dx p d 2 q' j dx K dx p dx p d 2 x' x dq' j \ 
\ dx' v dx" 7 dx p dx p dx w dx" 7 dx' x dx p dx K dx^ J 

Both formulas indicate that J 2 E and J 2 E are indeed affine bundles over E, with dif- 
ference vector bundles as stated above. 

The equivalence between the definition of the second order jet bundle given here 
and the traditional one is obtained observing that the iterated jet of a (local) 

section (p of E assume values not only in J 2 E but even in J 2 E, due to the Schwarz 
rule. Therefore, second order jets in the traditional sense, that is, classes of (local) 
sections where the equivalence relation is the equality between the Taylor expansion 
up to second order, are in one-to-one correspondence with these iterated jets of (local) 
sections. Moreover, a global section (p of E over M naturally induces a global section 
j 2 (p of J 2 E over M such that in the local coordinates used above 

fip{x) = (x p 1 ip i (x),d^ i (x),d fl d u ^ i (x)) . 

where d^ = d/dx p ] this is symbolically summarized by writing j 2 {p = (ip,dtp,d 2 (p). 



2.5 The Legendre Transformation 

A Lagrangian field theory is defined by its configuration bundle E over M and its 
Lagrangian density or simply Lagrangian, which in the present first order formalism is 
a map of fiber bundles over E: 

L-.^E — > f\ n T*M . (35) 
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The requirement that £ should take values in the volume forms rather than the func- 
tions on space-time is imposed to guarantee that the action functional S : T(E) — > R 
given by 

S[<f] = [ &((f,d<p) for V eT(E) (36) 
J At 

be well-defined and independent of the choice of additional structures, such as a space- 
time metric. 3 Such a Lagrangian gives rise to a Legendre transformation, which comes 
in two variants: as a map 

F£ : J X E — ► J l ®E (37) 

or as a map 

F£j : J X E — > J l ®E (38) 

of fiber bundles over E. For any point 7 in J\E, the latter is defined as the usual fiber 
derivative of £ at 7, which is the linear map from J\ E to A" T*M given by 

F£( 7 )-k = -^-£( 7 + A£) for k e J\E , (39) 

UiA A=0 

whereas the former encodes the entire Taylor expansion, up to first order, of £ around 7 
along the fibers, which is the affine map from J\E to A" T*M given by 

EC(7)-k = £(7) + - 7T £(7 + A(k-7)) for k^J\E . (40) 

d A A=0 

Of course, F£ is just the linear part of FL, that is, its composition with the bundle 
projection rj from extended to ordinary multiphase space: F£ = 77 o F£. In local 
coordinates as before, F£ is given by 

*-$"- £ -S< (4i) 

where L = L d n x. Finally, if £ is supposed to be hyperregular, which by definition 
means that F£ should be a global diffeomorphism, then one can define the De Donder- 
Weyl Hamiltonian "K to be the section of J 1@ E over J l ®E given by 

•K = F£ o (FL)- 1 . (42) 

In local coordinates as before, this leads to 

H = PX ~ L ( 43 ) 

where L = L d n x and "K = — H d n x, as stipulated in eq. (J3J). 

Conversely, the covariant Hamiltonian formulation of a field theory that can be 
described in terms of a configuration bundle E over M is defined by its Hamiltonian 

3 Strictly speaking, the integration should be restricted to compact subsets of space-time, which 
leads to an entire family of action functionals. 
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density or simply Hamiltonian, in the spirit of De Donder and Weyl, which in global 
terms is a section of extended multiphase space J 1 ®^ as an affine line bundle over 
ordinary multiphase space J 1 ®E: 

"K : J l ®E — ► J^E . (44) 

Such a Hamiltonian gives rise to an inverse Legendre transformation, which is a map 

FJ{ : J l ®E — ► J X E (45) 

of fiber bundles over E defined as follows. For any point z in J^®E, the usual fiber 
derivative of 0~C at z is a linear map from Jl®E to J\®E which when composed with 
the projection r\ from 3\®E to J]®E gives the identity on J]®E (since % is a section): 
such linear maps form an affine subspace of the vector space of all linear maps from 
Jl®E to Jl®E that can be naturally identified with the original affine space J\E, as 
explained in the Appendix. In local coordinates as before, FJi is given by 

4 - g <«> 

where "K = — H d n x. Finally, if "K is supposed to be hyperregular, which by definition 
means that F3i should be a global diffeomorphism, then one can define the Lagrangian 
L to be given by 

£(7) = (JCo(FM)- 1 )( 7 )-7. (47) 
In local coordinates as before, this leads to 

L = pH ~ H («) 

where L = L d n x and "K = — H d n x. 

Thus in the hyperregular case, the two processes are inverse to each other and 
allow one to pass freely between the Lagrangian and the Hamiltonian formulation. 
Of course, this is no longer true for field theories with local symmetries, in particular 
gauge theories, which require additional conceptual input. 

At any rate, it has become apparent that even in the regular case, the full power 
of the multiphase space approach to geometric field theory can only be explored if one 
uses the ordinary and extended multiphase spaces in conjunction. 



2.6 Canonical Forms 

The distinguished role played by the extended multiphase space is due to the fact that 
it carries a naturally defined multisymplectic form u, derived from an equally naturally 
defined multicanonical form 6 by exterior differentiation: it is this property that turns 
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it into the field theoretical analogue of the cotangent bundle of mechanics. 4 Global 
constructions are given in the literature [13-15], so we shall content ourselves with 
stating that in local coordinates q l ,p^ ,p) as before, 9 takes the form 

6 = pi dq l Ad n x tl + pd n x , (49) 

so uj = — d9 becomes 

uj = dq 1 a dp% a d n x ii — dp a d n x . (50) 

Given a Lagrangian &, we can use the associated Legendre transformation F£ to pull 
back 9 and uj and thus define the Poincare-Cartan forms 9 L and uj l on J X E associated 
with the Lagrangian L: 

C = (¥L)*9 , uj c = (¥L)*uj . (51) 

Similarly, given a Hamiltonian !K, we can use it to pull back 9 and uj and thus define 
the De Bonder- Weyl forms 9^ and uj^ on J 1 ®E associated with the Hamiltonian "K: 

M = <K*0 , uj m = 0K*uj . (52) 

Of course, uj l = —d6 L and uj^ = — d9< K ; moreover, supposing that "K o F£ = F£, 
we have 

e L = (wnyex , u L = (f£)Vk. (53) 

In local coordinates as before, eq. PH]) implies that 

BT / BT \ 

9 L = —dq^d n x,+ [L-—q;)d n x, (54) 



6^ = p f *dq i Ad n x^ - H d n x . (55) 

It is useful to note that the forms 9 L and 9^ allow us to give a very simple definition of 
the action functional: it is given by pull-back and integration over space-time. 4 Thus 
in the Lagrangian framework, the action associated with a section ip of E over M is 
obtained by taking the pull-back of 9 L with its derivative which is a section (ip, dip) 
of J X E over M, 

S[ip] = f (if,dip)*9 c for V cT(E) , (56) 

J M 

whereas in the Hamiltonian framework, the action associated with a section (ip, ir) 
of J 1 ®E over M is simply 

%>,tt] = h<p,*)*e K for {<p,n)eT{J 19 E) . (57) 

J M 

4 Note that this statement fails if one uses the ordinary duals instead of the twisted ones. 
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In both cases, it can be shown that the stationary points of the action are precisely 
the solutions of the corresponding Euler - Lagrange and De Donder-Weyl equations, 
respectively. It is therefore no surprise that these equations can be formulated globally 
through the vanishing of certain (in general nonlinear) differential operators £ and D 
defined solely in terms of the forms lj l and u^, respectively. However, an explicit 
construction in the spirit of global analysis [24,25] does not seem to be readily available, 
although there do exist various attempts that go a long way in the right direction; see, 
e.g., [26] for the Lagrangian case and [13] for the Hamiltonian case. 



2.7 Euler - Lagrange and De Donder-Weyl Operator 

Theorem 1 Given a Lagrangian density as in eq. 4,9,5j) above, define the corresponding 
Euler - Lagrange operator to be the map 

l:.J 2 E — > V*E®/\ n T*M (58) 

of fiber bundles over J X E 5 that associates to each 2-jet (p,dp,d 2 p) of (local) sections 
ip of E over M and each vertical vector field V on E the n-form on M given by 

E(<p, dp, d 2 cp) - V = (<p, dp)* (i Jly u c ) . (59) 

Then for any section p> of E, E,((p,d(p,d 2 p) is the zero section if only if (p satisfies the 
Euler - Lagrange equations associated to £ . 



PROOF: Let V be a vertical vector field on E, with local coordinate expression 

• d 

v = v 'w 

(cf. eq. (fT2J)). and let J X V be its lifting to J X E, with local coordinate expression 



Pint \ Pink ^ 



dq 1 \ dq k ^ dx^ J cM 

(cf. eq. (|T3|) ). Applying the exterior derivative to eq. (jSU), contracting with J X V and 
then pulling back with (ip, dp) gives, after some calculation, 

(p,dpy(i Jlv u L ) = ^-(p,dp)v\p)d n x + ^-(p,dp)v\p)dyd n x 

d 2 T 8T 
+ 1 r- n -(^dp)V\p)d^ d n x - — {p,dp)V\p)d n x 
oqi dq 1 ^ oq 



3 Again, we suppress the symbols indicating the pull-back of bundles from E or M to J 1 .©. 
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where it is to be noted that the terms depending on the derivatives of V have dropped 
out. This leads to the following explicit formula for £: 

£&,dp,d\) = ^(ZLfadtpj) - f^, W) <tf®d n x . (60) 

In particular, it is clear that £ depends on p> only through the point values of ip and 
its partial derivatives up to second order, which concludes the proof. □ 



Theorem 2 Given a Hamiltonian density as in eq. \44\) above, define the corresponding 
De Donder- Weyl operator to be the map 

V:J\,P®E) — ► V*(.P®E)® /\ n T*M (61) 

of fiber bundles over J l ®E 6 that associates to each 1-jet (ip, n, dip, dir) of (local) sections 
(p, 7r) of J 1 ®E over M and each vertical vector field V on J l ®E the n-form on M given 
by 

2%, 7T, dp, dn)-V = (p, vr)* (i y Wm) • (62) 

Then for any section ((p, n) of J 1 ®E, T)(p, n, dp, dn) is the zero section if only if (<p, n) 
satisfies the De Donder- Weyl equations associated to "K. 

PROOF: Let V be a vertical vector field on J 1 ®E, with local coordinate expression 

• d d 
V = V 1 -7— + V? 



dq i ' 1 dpf 

Applying the exterior derivative to eq. ()55|) . contracting with V and then pulling back 
with (<p, 7r) gives, after a short calculation, 

dH 

(p,KY(i v ujk) = dXV r (p,n)d n x + —(<p,ir)V i (<p,ir)d n x 



dq l 

dH 

dp- 



d^ Vf(<p,ir) d n x + ^Oatt) Vf(<p,ir) d n x . 



This leads to the following explicit formula for D: 

V(<p,ir,d<p,dir) = (J^,70 - drf^ dq* ® d n x 



dH 

+ [ Qjfifa'K) + <W ) dp^d n x 



(63) 



In particular, it is clear that D depends on (p, it) only through the point values of p 
and 7r and their partial derivatives up to first order, which concludes the proof. □ 

6 Again, we suppress the symbols indicating the pull-back of bundles from E or M to J 1 ®E. 
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2.8 Jacobi Operators 



In order to make contact with the functional formalism to be discussed in the next 
section, we must also derive explicit expressions for the linearization of the Euler- 
Lagrange operator and the De Donder-Weyl operator around a given solution of the 
equations of motion. This leads to linear differential operators between vector bundles 
over M that we shall refer to as Jacobi operators, generalizing the familiar derivation 
of the Jacobi equation by linearizing the geodesic equation. 

In its Lagrangian version, the Jacobi operator is a second order differential operator 

3M ■ r(vy — r(v®) , (64) 

where V v = <p*{VE) and V® = tp*{V*E) <g> f\ n T*M, obtained by linearizing the 
Euler - Lagrange operator £ around a given solution tp of the equations of motion. 
Similarly, in its Hamiltonian version, the Jacobi operator is a first order differential 
operator 

a^^rxw — > r(v® n) ), (65) 

where V (v , v) = (tp, tt)*(V( J 1 ®^)) and V ( ® w) = (tp, 7r)*(V* (J 1 ®E))®/\ n T*M, obtained 
by linearizing the De Donder-Weyl operator T> around a given solution (tp,n) of the 
equations of motion. (Thus in both cases, the vector bundles involved are obtained by 
pulling back the appropriate vertical bundle and its twisted dual with the solution of 
the nonlinear equation around which the linearization is performed,) To obtain explicit 
expressions, consider an arbitrary variation tp x around tp and evaluate £(ip\, dtp\, d 2 tp x ) 
which, for each A, is a section of V® x , observing that since tp = tp\\ x=Q is a solution, 
£(tp\, dcpx, d 2 tp x )\ x _ is the zero section of V®, and setting 



A 9 

**P = ^A 



A=0 



(66) 



Noting that in local coordinates, the value of £(<^a, dtp\, d 2 tp x ) at a point x in M with 
coordinates x M has coordinates (x^, tp\(x), 8,(tp x , dtp x , d 2 tp x )i(x)) where the last piece is 
the coefficient of dq % <E> d n x in eq. (JBDJ), we get by differentiation with respect to A 



d 

— £(tp x ,dtp x ,d 2 tp x ) 
x i 9 

= 



A=0 



d 2 L d 2 L \ 

-(if, dtp) b~tp> : {tp, dtp) d v btp? ) dq l ® d n x . 
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Similarly, consider an arbitrary variation (<p\, tt\) around (ip, tt) and evaluate 
D(ip\, 7T\, d(f\, dwx) which, for each A, is a section of V, 



((p,n) = (px,tt x ) 
and setting 



© 



=0 is a solution, D(<fx,^x,dipx,dTT X ] 



observing that since 
I . . is the zero section of V® n , 

lA=0 (vw' 



d 

(6<p,8ir) = —(<fx,^x) 



A=0 



(67) 



Again, noting that in local coordinates, the value of D((p\, 7t\, d<p\, dux) at a point x 
in M with coordinates x M has coordinates (p\(x), (ir\)i(x), D({p\, tt\, d(f\, dit\)i(x), 
D(<p\, 7Ta, dipx, dnxY^x)) where the last two pieces are the coefficients of dq 1 ®d n x and 
of dp? ® d n x in eq. ((H3J), we get by differentiation with respect to A 

d 

X>(^A,7r A ,<9y3A,<97r A ) 

A=0 

a 



+ 



+ 



d 2 H 
dqi dq i 

d 2 H 
dqi dpi 



-0,vr) 5ip> + 



d 2 H 

dp V j dq 

d 2 H 

dp" dp\ 



{<f,7T) 8ir) 



dq 1 <g> d n x 
(ip,7r) 5< + 9„ V) dpt®d n x . 



In order to show how to extract the Jacobi operators from these expressions, by 
means of a globally defined prescription, we apply the following construction [10]. 
Let F be a fiber bundle over M, with bundle projection ttf,m '■ F — > M, and W 
be a vector bundle over F with bundle projection ~k-w,f '■ W — > F, which is then also a 
fiber bundle (but not necessarily a vector bundle) over M with respect to the composite 
bundle projection nw,M = ^f,m ° ^w,f '■ W — > M. Thus W admits two different kinds 
of vertical bundles, VpW and VmW, with fibers defined by (Vf) w W = ker T w ttw,f and 
(Vm)wW = ker T w n w>M for w e IV; obviously, the former is contained in the latter as a 
vector subbundle. Moreover, since W is supposed to be a vector bundle over F, there 
is a canonical isomorphism VfW = n^pW. On the other hand, consider the vertical 
bundle VF of F which can be pulled back to W to obtain a vector bundle 7r^ F (VF) 
over W, with fibers defined by (tt^ f (VF)) w = VfF = ker Tfixp,M for weW with 
/ — n w,FW. Note also that the tangent map to the bundle projection tiw,f, which 
by definition has kernel VpW, maps VmW onto VF, so we have the following exact 
sequence of vector bundles over W: 

— > VpW tt^ f IV — > VmW — 7r^ F (VF) — > . 

The crucial observation is now that this exact sequence admits a canonical splitting 
over the zero section : F — > IV, given simply by its tangent map. Indeed, its 
tangent map TfO : TfF — ► T (j)IV a ^ an y point /eF takes the vertical subspace 
V/F to the M-vertical subspace (Vm)o(/)W / and so restricts to a vertical tangent map 
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VfO : VfF — > (Vm)o(/)W / whose composition with the restriction of the tangent map 
Tq(j)Hw,f '■ To(f)W — > T/F to (Vm)o(/)W gives the identity on V/F. Thus the image of 
VfO is a subspace of (Vm)o(/)W that is complementary to the subspace (Vf)o(/)W^ — Wf 
and provides a surjective linear map 0$ : (Vm)o(/)W / - ► W/ °^ which it is the kernel. 
At the level of bundles, this corresponds to a surjective vector bundle homomorphism 
o : V M W\ Q - W. 

Applying this construction to the situation at hand, take F = E in the Lagrangian 
case and F = J 1@ E in the Hamiltonian case, setting W = V*(F) [\ n T*M in 
both cases. The fact that the Euler - Lagrange or De Donder-Weyl operator is being 
linearized around a solution tp or (</?, it) of the equations of motion then means that 
we are evaluating its derivative, which a priori takes the variation Sip or (5<p, 5n) to a 
vector field on W along M which is vertical with respect to the projection of W onto M, 
precisely over the zero section, so we can apply the operator a just introduced to project 
it down to a section of W over M itself. This operation completes the definition of the 
Jacobi operators, namely 

= cr(^Z{ Vx ,dp x ,d 2 p x ) ) , (68) 



A=0 



(69) 



ana 

3 M [(p, 7r] • (Sip, 5n) = a (J^T>((p x ,Trx,d<px,dir x ) ^ . 
The local coordinate expressions are the ones derived above, that is, 

dM ■ S( P = { aJ a (y» d( P) d^djipj 
I oqt dq 1 ^ 

+ {<ik ( ^ + & ~ ww>' dv) ) d ^ (70) 

and 

= {wk^ v + *** " d > 6< ) dql dn% (71) 



20 



3 Functional Approach 



Let us begin by recalling the definition of the Poisson bracket between functions on a 
symplectic manifold with symplectic form u. First, one associates to each (smooth) 
function / a (smooth) Hamiltonian vector field Xf, uniquely determined by the condi- 
tion 

ix f u = df . (72) 

Then the Poisson bracket of two functions / and g is defined to be the function {/,<?} 
given by 

{f,g} = -i Xf i Xg u; = df(X g ) = -dg(X f ). (73) 

The goal of this section is to show that formally, the same construction applied to 
covariant phase space links the Witten symplectic form to the Peierls bracket. 

3.1 Covariant Phase Space 

In contrast to the traditional non-covariant Hamiltonian formalism of field theory, where 
phase space is a "space" of Cauchy data, covariant phase space, denoted here by §, is the 
"space" of solutions of the equations of motion, or field equations. Of course, one cannot 
expect these two interpretations of phase space to be equivalent in complete generality, 
since it is well known that, for nonlinear equations, time evolution of regular Cauchy 
data may lead to solutions that, within finite time, develop some kind of singularity. 
An even more elementary prerequisite is that the underlying space-time manifold M 
must admit at least some Cauchy surface £ : this means that M should be globally 
hyperbolic. 

Thus our basic assumption for the remainder of this paper will be that the underlying 
space-time manifold M should be globally hyperbolic. Globally hyperbolic space-times 
are the natural arena for the mathematical theory of hyperbolic (systems of) partial 
differential equations, in which the Cauchy problem is well posed. There are by now 
various and apparently quite different definitions of the concept of a globally hyperbolic 
space-time, but they have ultimately turned out to be all equivalent; see Chapter 8 
of [27] for an extensive discussion. For our purposes, the most convenient one is that M 
admits a global time function whose level surfaces provide a foliation of M into Cauchy 
surfaces, providing a global diffeomorphism M = R x E. As an immediate corollary, 
we can define the concept of a (closed/open) time slice in M: it is a (closed/open) 
subset of M which under such a global diffeomorphism corresponds to a subset of the 
form I x E where I is a (closed/open) interval in R. 

In the Lagrangian as well as in the Hamiltonian approach to field theory, the equa- 
tions of motion are derived from a variational principle, that is, their solutions are the 
stationary points of a certain functional S called the action and defined on a space 
of sections of an appropriate fiber bundle over space-time which is usually referred to 
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as the space of field configurations of the theory and will in what follows be denoted 
by C. More concretely, C is the space T(F) of smooth sections of a fiber bundle F 
over M: in the Lagrangian approach, F is the configuration bundle E, whereas in the 
Hamiltonian approach, F is the multiphase space J 1 ®E, regarded as a fiber bundle 
over M. 

Formally, we shall as usual think of C as being a manifold (which is of course infinite- 
dimensional). As such, it has at each of its points a tangent space T^G that can be 
defined formally as a space of smooth sections, with appropriate support properties, of 
the vector bundle = <p*(VF) over M, i.e., T^Q c r°°(V^). The cotangent space T^Q 
will then be the space of distributional sections, with dual support properties, of the 
vector bundle V® = <p*(V*F) f\ n T*M over M, i.e., T^C c Y-°°(V®). It contains 
as a subspace the corresponding space of smooth sections, where the pairing between 
a smooth section of and a smooth section of V? (with appropriate support condi- 
tions) is given by contraction and integration of the resulting form over M. Similarly, 
the second tensor power T^C ® T^C of T^C can be thought of as the space of distribu- 
tional sections, again with dual support properties, of the second exterior tensor power 7 
V® M V® of V®; it contains as a subspace the corresponding space of smooth sections, 
where the pairing between a pair of smooth sections of and a smooth section of 
V®MV® (with appropriate support conditions) is given by contraction and integration 
of the resulting form over M x M. 

Regarding the support conditions to be imposed, the first two options that come 
to mind would be to require that either the elements of T^C or the elements of T^C 
should have compact support, which would imply that the support of the elements of 
the corresponding dual, T^C or T^C, could be left completely arbitrary: 

Option 1 : Tfi = r°°(V*) , T;Q = T;°°(V®) (74) 

Option 2 : T^C = r~(V*) , T;Q = T-°°(V®) (75) 

There is a third option that makes use of the assumption that M is globally hyperbolic. 
To formulate it, we introduce the following terminology. A section of a vector bundle 
over M is said to have spatially compact support if the intersection between its support 
and any (closed) time slice in M is compact, and it is said to have temporally compact 
support if its support is contained in some time slice. Then, as in Ref. [8], we require 
the elements of T^C to have spatially compact support and the elements of X^C to have 
temporally compact support: 

Option 3 : 7>C = r-(^) , T;e = T^(V®) (76) 

Obviously, for each of these three options, the two spaces listed above are naturally 
dual to each other. 8 

7 If V and W are vector bundles over M, V 23 W is denned to be the vector bundle over M x M 
with fibers given by (VH W) {x ^ v) =V x ®W y , for all x, y e M. 

8 Here and in what follows, the symbols T c , T sc and Y tc indicate spaces of sections of compact, 
spatially compact and temporally compact support, respectively. 
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These constructions can be applied to elucidate the nature of functional derivatives 
of functionals on C, such as the action. Namely, given a (formally smooth) functional 
F : C — > R, its functional derivative at a point is the linear functional on T^C which, 
when applied to 50 , yields the directional derivative of Fat along 50 , defined by the 
requirement that for any one-parameter family of sections 0a of F such that <p\ \ 



I A=0 



F'[4>] ■ 50 = ^ F[0 A ] 



if 



A=0 



A=0 



Then F'[0] is a distributional section of with appropriate support properties (dual 
to those required for T^C). In local coordinates, its action on 50 can (formally and at 
least when the intersection of the two supports is contained in the coordinate system 
domain) be written in the form 



F'[0]-50 = f d n x %[<t>]{x) ■ 5<f>(x) , 



(77) 



The expression (5F/50)[0], sometimes called the variational derivative of Fat , is then 
a distributional section of V£ (over the coordinate system domain). In the Lagrangian 
framework, 

whereas in the Hamiltonian framework, 



5F 

50 



Similarly, the second functional derivative of Fat is the symmetric bilinear functional 



on T^C which, when applied to 5(f) l and 
for any two-parameter family of sections 



)0 2 , can be defined by the requirement that 
b XlM of F such that ^Ai,a 2 | A1)A2=0 = 0, 



F"\ 



d 2 



d\ 1 d\ 2 



3 Ai,A 2 J 



Ai,A 2 =0 



if 



Hi 



d 



0Ai, 



A 2 



Ai,A 2 =0 



50 2 



_d_ 

OX, 



0Ai 



A 2 



Ai,A 2 =0 



Then F"[0] is a distributional section of V®MV® with appropriate support properties 
(dual to those required for T^C ® T^C). In local coordinates for M x M induced from 
local coordinates for M, its action on (<50 1 ,50 2 ) can (formally and at least when the 
intersection of the supports is contained in the coordinate system domain) be written 
in the form 



F"[0]- (50^) 



d " x 1 ^ T&W( x >y) ■ WiWMziy)) , (78) 
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The expression (5 2 F/5(fi 2 )[(fi], sometimes called the variational Hessian of F at <fi , is 
then a distributional section of VI Kl V£ (over the coordinate system domain). In the 
Lagrangian framework, 

5 2 F 5 2 F 

w mx ' v) = ww [ip]{x,y)dql ® dqJ 7 

whereas in the Hamiltonian framework, 



S 2 F 5 2 F S 2 F 

J^W&y) = Sipifiipj ^rffay) dq l ®dq J + j~qr^ [<P, A \?,V) dq % ® dp) 

5 2 F , 5 2 F 



+ ^ S(fij [<P, A {x, y) dp? g> dq 3 + — [cp, ir]{x, y) dp? ® dp) 



Of course, for the integrals in eqs (J77)) and (J78j) to make sense, even when interpreted 
in the sense of pairing distributions with test functions, we must make some assumption 
about support properties, which leads us back to the options stated in eqs (|71 |) -(|7H |) . 
Option 1: when F is arbitrary, we have to restrict the sections 5cfi, 8(f) 1 , 5<fi 2 of V<j> 
considered above to have compact support (which can be achieved if the sections (fix, 
4>x u \ 2 of F are supposed to be independent of the parameters outside a compact subset). 
Option 2: when F is local, which we understand to mean that its functional dependence 
on the fields is non-trivial only within a compact region, or equivalently, that its func- 
tional derivative F'[<f>] at each (fi has compact support, the sections 5 (ft, 5(fi 1 , 5<fi 2 of 
considered above may be allowed to have arbitrary support; this is the case for local 
observables defined as integrals of local densities over compact regions of space-time 
and, in particular, over compact regions within a Cauchy surface E (energy, momentum, 
angular momentum, charges etc. within a finite volume). Option 3: when F is local 
in time, which we understand to mean that its functional dependence on the fields is 
non-trivial only within a time slice, or equivalently, that its functional derivative F'[<f>] 
at each <fi has temporally compact support, we have to restrict the sections 5(f>, S(f> 1 , 
5<p 2 of considered above to have spatially compact support (which can be achieved 
if the sections (fix, <px 1 ,x 2 are supposed to be independent of the parameters outside a 
spatially compact subset); this is the case for global observables defined as integrals of 
local densities over time slices and, in particular, over a Cauchy surface S (total energy, 
total momentum, total angular momentum, total charges etc.). 

Finally, covariant phase space S is defined to be the subset of C consisting of the 
critical points of the action: 

S = {<fiee/S'[<fi} = 0} . (79) 



Formally, we can think of S as being a submanifold of C whose tangent space at any 
point of § will be the subspace T^S of the tangent space T^C consisting of the solu- 
tions of the linearized equations of motion (where "linearized" means "linearized around 
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the solution of the full equations of motion"), which are precisely the sections of 
belonging to the kernel of the corresponding Jacobi operator S[4>] '■ r(V^) — > T(V®) : 

T^S = ker 3[4>] ■ (80) 



3.2 Symplectic Structure 

Our next goal is to justify the term "covariant phase space" attributed to S by show- 
ing that, formally, S carries a naturally defined symplectic form Q, derived from an 
equally naturally defined canonical form by formal exterior differentiation. Accord- 
ing to Crnkovic, Witten and Zuckerman [1-3] (see also [4]), the symplectic form Q 
can be obtained by integration of a "symplectic current", which is a closed (n — 1)- 
form on space-time, over an arbitrary spacelike hypersurface S. Here, we show that 
this "symplectic current" can be derived directly from the multisymplectic form u or, 
more explicitly, from the Poincare- Cartan form ujh in the Lagrangian approach and 
the De Donder- Weyl form ujt in the Hamiltonian approach. 

We begin with the definition of B and fl in terms of 9 and u, which is achieved by 
a mixture of contraction and pull-back: given a point in C (a smooth section of F) 
and smooth sections 5(f), S(j) l , 5<f> 2 of V^, insert 5(f) into the first of the n arguments of 9 
or 54>\ and 5(f) 2 into the first two of the n + 1 arguments of u and apply the definition 
of the pull-back with (which amounts to composition with the derivatives <90 of 0) 
to the remaining n — 1 arguments to obtain (n — l)-forms on space-time which are 
integrated over S. Note that these integrals exist if we assume that 5(f) and either 5<p 1 
or 5(f> 2 have spatially compact support, since this will intersect S in a compact subset. 

Explicitly, in the Lagrangian framework, we have 



6^(50) = JtoWOtfadSv) (81) 

and 

fi (50i,<50 2 ) = J (<p,d<p)* u c (6ipi,d 8(fx , 5^ 2 ,d5f 2 ) (82) 

where the notation is the same as that employed in eq. (JoT)j) : = (p is a section of 
E over M and j 1 ^ = (if , dtp) is its (first) prolongation or derivative, a section of J l E 
over M, while 5<f> = Sip, 5<f> 1 = 8ip l , 5<f> 2 = 5ip 2 are variations of = ip, all sections 
of VE over M, and <5jV = (8ip,dSip), Sj 1 ip 1 = (Sip^d Sip^, 5j 1 ip 2 = (Sip 2 ,d Sip 2 ) are 
the induced variations of j x ip = (ip,dip), all sections of V(J 1 E) = ^(VE) over M. 
In local coordinates, 



A 9 

6<p = Tx f, 



a=o oq % 



and 



.•9 ~ . , d 



5ip l — + d^Stp* 



a=o dq 1 dq 



A' 
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whereas 6 L is given by eq. and uj l by 



10, 



+ 



Qqj Q g i ^ A ^ + o,.- ^ A ^ ) A rf% 

<9 2 L 



cte^ 1 dq^ 



dql dql 
dq % a d n x — d(L - 



dL 



dq l ^ 



q % ) a (i n x . 



(The exterior derivative in the last term could be worked out explicitly, but we shall 
not need this expression because the last two terms vanish under contraction with two 
vertical vectors.) Then 



a 



and 



with the "symplectic current" J given by 



1? "^2; 



dqi dq 



(<P, df) (5<P\ <V 2 - V2 Vi) 



+ 



A' 

<9 2 L 



or equivalent ly 



i5 u r2J 



d 2 L d 2 L 

(p,9ip) Vi + ttt" O, dtp) d v 6<p[ ) V 2 



9g J dq 



dql dq % 



13) 



H) 



15) 



(86) 



The same results can be obtained even more directly in the Hamiltonian framework, 
in which we have 



and 



15 u r2) 



(<p, 7r)* 9h{Sp, Sir) 



(ip, vr)* u M {Spi, fax , 5ip 2 , 5tt 2 ) 



i7) 



where the notation is the same as that employed in eq. (JoTj) : <fi = (ip, it) is a section of 
J l ®E over M while 50 = (5<p, Sn), 5<p 1 = (Sip 1 , Sn-^), <50 2 — (^25 are variations of 
cf) = (p,7r), all sections of V(J 1 ®E) over M. In local coordinates, 



^ = ex** 



A=0 



^ a? ' 57r = aT A 



<9 



A=0 
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whereas 9^ is given by eq. (J55|) and by 

u)<}t = dq l a dpi a d n x^ — dH a d n x 

(The exterior derivative in the last term could be worked out explicitly, but we shall not 
need this expression because the last term vanishes under contraction with two vertical 
vectors.) Then 



6^(50) = / da M vrf V (89) 

and 



E 



with the "symplectic current" J given by 

4(5^,5^) = 6<p\ Sirrf - 6<pi Sir J . (91) 

Incidentally, these formulas show that, just like in mechanics, the canonical form O and 
the symplectic form f2 do not depend on the choice of the Hamiltonian Ji. 

Another important result, duly emphasized in the literature [1-4], is the fact that 
on covariant phase space S, the symplectic form Q does not depend on the choice of the 
hypersurface S used in its definition, since for any solution <p of the equations of motion 
and any two solutions 5<fi 1 , <5</> 2 of the linearized equations of motion, the "symplectic 
current" J ( f > (5(p 1 , 50 2 ) is a closed form on space-time. To prove this, assume that is a 
point in § and observe that a tangent vector 5(f) in T^C belongs to the subspace T^S if 
and only if 5(j), as a section of V&, satisfies the pertinent Jacobi equation, which reads 



/ d 2 L d 2 L \ 

in the Lagrangian framework and 

d 2 H d 2 H 



(92) 



(93) 



in the Hamiltonian framework. Thus if S(f) 1 and <50 2 both satisfy the Jacobi equation, 
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we have 



+ 



+ (^ (v? ' 9v?) V2 + ^% ( ^ } w Wi 

in the Lagrangian framework and 

+ + af^'*' fa « ) 

- ^ ( eff^' *> S ' A + £fk> (v ' n) s *» ) 

in the Hamiltonian framework: obviously, both of these expressions vanish. 

Of course, independence of the choice of hypersurface holds only for but not for 0. 
In fact, if M lj2 is a region of space-time whose boundary is the disjoint union of two 
hypersurfaces Si and E 2 , then Q^, 2 = but 

e S2 -e Sl = 5s Ml , 2 (94) 

where Sm 1i2 is the action calculated by integration over M 12 and 5 is the functional 
exterior derivative, or variational derivative, on S. 
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3.3 Poisson Bracket 



Given a relativistic field theory with a regular first-order Lagrangian, one expects each 
of the corresponding Jacobi operators S[4>] (0 e §) to form a hyperbolic system of second- 
order partial differential operators. A typical example is provided by the sigma model, 
where E is a trivial product bundle M x Q, with a given Lorentzian metric g on the 
base manifold M, as usual, and a given Riemannian metric h on the typical fiber Q. 
Its Lagrangian is 

so that the coefficients of the highest degree terms of the Jacobi operator 3[<p] are 



dql dq^ 

which clearly exhibits the hyperbolic nature of the resulting linearized field equations. 

A general feature of hyperbolic systems of linear partial differential equations is the 
possibility to guarantee existence and uniqueness of various types of Green functions. 
In the present context, what we need is existence and uniqueness of the retarded Green 
function G^, the advanced Green function G^ and the causal Green function G^ for 
the Jacobi operator £f[0], for each <f> e §. By definition, the first two are solutions of the 
inhomogeneous Jacobi equations 

3MGi(x,y) = S(x,y) , 3 y [<!>] G$(x, y) = S(x,y) , (95) 

or more explicitly, 

2M km G± m \x,y) = 5 l k 5(x,y) , 3y[^\k m Gf m {x,y) = 5 l k 5(x,y) , (96) 

where dz[4>] denotes the Jacobi operator with respect to the variable z, characterized by 
the following support condition: for any x,yeM, G^(x,y) = when x£j + {y) and 
G^(x,y) = when x<£ J~(y), where J + (y) and J~(y) are the future cone and the past 
cone of y, respectively. The causal Green function, also called the propagator, is then 
simply their difference: 

G^ = G^ - G+ . (97) 
Obviously, it satisfies the homogeneous Jacobi equations 

dxWG+fay) = , dyWG+fay) = 0. (98) 

Note that the symmetry of the Jacobi operator d[(f>], stemming from the fact that it 
represents the second variational derivative of the action, forces these Green functions 
to satisfy the following exchange and symmetry properties: 

Gf\y )X ) = Gf\x,y) , G l *(y,x) = -G k J(x,y). (99) 
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It should be pointed out that existence and uniqueness of these Green functions cannot 
be guaranteed in complete generality: this requires not only that M be globally hyper- 
bolic but also that the linearized field equations should form a hyperbolic system. Here, 
we shall simply assume this to be the case and proceed from there; further comments 
on the question will be deferred to the end of the section. 

Our next step will be to study certain specific (distributional) solutions X F [(f>] of the 
general inhomogeneous Jacobi equation 

= F'[<f>] (ioo) 

for smooth functionals F on covariant phase space which are (at least) local in time. 
To eliminate the ambiguity in this equation stemming from the fact that the functional 
derivative F [(f)] on its rhs belongs to the space T^S which is only a quotient space of 
the image space T^C of the Jacobi operator S[4>] (an inclusion of the form T^S cT^C 
induces a natural projection from T^C to TTS), it is necessary to first of all extend 
the given functional F on § to a functional F on C of the same type (smooth and 
local in time), whose functional derivative F'[0] does belong to the space T^C which, 
as we recall from eq. 1)760. consists of the distributional sections of V® of temporally 
compact support. Next, convolution with the retarded and advanced Green function 
introduced above produces formal vector fields over § which to each solution cf) e § of the 
field equations associate (distributional) sections X~^[<p] and X^[</>] of V^, respectively. 
In local coordinates, their definition can (formally and at least when the intersection 
of the two supports is contained in the coordinate domain) be written in the form 

X%<P]\x) = f d B yGj"(a;,y)^M(y). (101) 

Both of them satisfy the inhomogeneous Jacobi equation 

SmX$[<P}) = F'[0] . (102) 

Similarly, convolution with the causal Green function leads to another formal vector 
field over § which to each solution e S of the field equations associates a (distribu- 
tional) section X^[<f>] of V^. Again, in local coordinates, its definition can (formally 
and at least when the intersection of the two supports is contained in the coordinate 
domain) be written in the form 

= / AG^J/jS^). (103) 

It satisfies the homogeneous Jacobi equation 

0[0](X^>]) = 0, (104) 

since according to eq. ()97)l 

X&\ = X^]-X±[0] . (105) 
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Note that the convolutions in eqs (jlUlj) and ()103p exist due to our support assumptions 
on F (requiring F'[(f>] to have temporally compact support) and due to the support 
properties of the Green functions G^ and G^. 

According to eq. ()104|) . the prescription of associating to each solution 0eS of the 
field equations the section Xb[4>] of defines a formal vector field on S which is tangent 
to S . (It becomes more than just a formal vector field if Fis such that X^[(f)] belongs to 
T^S , which requires it to be not just a distributional section but a smooth section of 
and to satisfy appropriate support properties; we shall come back to this point later on.) 
The main statement about this formal vector field, to be proved below, is that (a) it 
does not depend on the choice of the extension F of F, so we may simply denote it 
by X F [(j)], and (b) that it is formally the Hamiltonian vector field associated to Fwith 
respect to the symplectic form Q discussed in the previous subsection. More explicitly, 
we claim that for any solution <fi e S of the field equations and any smooth section 5(f> 
of V,/, with spatially compact support which is a solution of the linearized field equations, 
we have 

^(X F [0],50) = F'[4>}-64>. (106) 

Note that under the assumptions stated, both sides of this equation make sense although 
we have originally defined Q ( / > (5(f> 1 , 5(f> 2 ) only in the case where both 5(f> 1 and 5(f> 2 are 
smooth; the extension of this definition, given in the previous subsection, to the case 
where one of them is a distribution is straightforward. 

To prove this key statement, let us begin by recalling that the symplectic form Q 
and the symplectic current J of the previous subsection are really defined on C and not 
only on S - the only difference is that on C, Q is only a presymplectic form so that J 
should be more appropriately called the presymplectic current and that J on C is no 
longer be conserved so that Q on C will depend on the choice of the hypersurface S. 
At any rate, we can almost literally repeat the calculation performed at the end of the 
previous subsection, either in the Lagrangian or in the Hamiltonian formulation, to 
show that for any solution e S of the field equations and any smooth section 5(f) of 
with spatially compact support, we have 

d^{x%<j>u<j>) = {Mu^W) W k - WW 1 ) x$i<p} k , (107) 

so that if 5(f) is a solution of the linearized field equations, 

fyJtiXpMM) = j^[<f>} ■ 5( t> ■ (108) 

Now since, by assumption, the support of (5F/S(j))[(j)] is contained in some time slice, 
we can choose two Cauchy surfaces £_ to the past and S + to the future of this time 
slice and, using that 5(f) has spatially compact support, integrate eq. (|108j) over the time 
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slice S- between £_ and S and similarly over the time slice S + between £ and S + . 
Applying Stokes' theorem, this gives 



fda^x) J£(Xj[<j>],6<j>)(x) = J da^x) Jf{X^[P]M)(x) + J d n x • 6<j>(x) , 

Ja,(x) J^X+^S^x) = jf da,(x) J$(X + ~[<p],5<f>)(x) + jf rf"a; ^[0](x) • 50(a;) . 

But the support conditions on G^, together with the fact that the support of (5F/ 5(p) [cp] 
lies to the future of £_ and to the past of £+, imply that X^[<f>\ vanishes on £_ and 
similarly that vanishes on £+, so the first term on the rhs of each of these 

equations is zero. Thus taking their difference and inserting eq. (|105jl . we get 



and since (6F/6<p)[<p] vanishes outside S^uS + , 

QtiXf&lSfi = [ d n x 5 -Lu]{ x ) ■ 5^){x) . (109) 
hi H 

Finally, observe that since 8<p is supposed to be a solution of the linearized field equa- 
tions (and hence tangent to §), the rhs of this equation does not depend on the choice 
of the extension F of F. Therefore, X^[<p] will not depend on this choice either provided 
the symplectic form is weakly non-degenerate. Now using the space-time split of M 
over £ provided by the tangent vector field d t of some global time function i on M 
or its dual dt, and identifying solutions 5<p of the linearized field equations with their 
Cauchy data on £, 9 it can be seen by direct inspection, either of eqs an d (jSSj) in 
the Lagrangian formalism or of eqs (fHU|) and (jHTjl in the Hamiltonian formalism, that 
the expression Q^{5(pi, Sfa) can only be zero for all S(p2 if 8<fi\ vanishes, as soon as we 
require the Lagrangian L to be regular in time derivatives, that is, to satisfy 

d 2 L 

det : r ^ , (110) 

dq dq J 

or equivalently, the Hamiltonian to be regular in timelike conjugate momenta, that is, 
to satisfy 

Moreover, it can be shown that this statement will remain true if 8(j)\ is allowed to be a 
distributional solution of the linearized field equations with arbitrary support, as long 



9 Explicitly, in the Lagrangian formalism, the Cauchy data for Sip on M are Sip and Sep on E, whereas 
in the Hamiltonian formalism, the Cauchy data for (Sip, Sir) on M are Sip and Sir on E. 
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as 5(j)2 runs through the space of smooth solutions of the linearized field equations with 
spatially compact support. 

Let us summarize this fundamental result in the form of a theorem. 

Theorem 3 With respect to the symplectic form Q on covariant phase space as defined 
by Crnkovic, Witten and Zuckerman, the Hamiltonian vector field X F associated with a 
functional F which is local in time is given by convolution of the functional derivative 
of F with the causal Green function of the corresponding Jacobi operator. 

Note that in view of the regularity conditions employed to arrive at this conclusion, 
the previous construction does not apply directly to degenerate systems such as gauge 
theories: these require a separate treatment. 

Having established eq. (jlOfij) . it is now easy to write down the Poisson bracket of 
two functionals Fand G on S : it is, in complete analogy with eq. (JZHJ), given by 

{F,G}[<P} = F'[0].X G [0] = -G>]-X F [0], (112) 

or 

{F , G}[0] = / d n x ^M(x) X G [<P]\x) = - / d n x ^ k [<P]{x) X F [<P]\x) . (113) 
JM 0( P JM 0( P 

Inserting eq. ()103|) . we arrive at the second main conclusion of this paper, which is an 
immediate consequence of the first. 

Theorem 4 The Poisson bracket associated with the symplectic form Q on covariant 
phase space as defined by Crnkovic, Witten and Zuckerman, according to the standard 
prescription of symplectic geometry, suitably adapted to the infinite- dimensional setting 
encountered in this context, is precisely the field theoretical bracket first proposed by 
Peierls and brought into a more geometric form by DeWitt: 

{F,G}[<P] = [ d n x [ d n y^}{x)G k ;{x,y)^[<P]{y). (114) 

JM JM 0( P 0< P 

Of course, for the expressions in eqs (|112|l - (|114jl to exist, it is not sufficient to require F 
and/or G to be local in time. In fact, if we want to use conditions that (a) are sufficient 
to guarantee existence of this Poisson bracket without making use of specific regularity 
and support properties of the propagator, (b) are the same for F and G and (c) are 
reproduced under the Poisson bracket, we are forced to impose quite rigid assumptions: 
the functionals under consideration must be assumed to be both regular and local, in 
the sense that their functional derivative at any point <p of S must be a smooth section 
of V® of compact support (this will force the corresponding Hamiltonian vector field 
to be a smooth section of V$ of spatially compact support). 



33 



On the other hand, it must be pointed out that this Poisson bracket, which we might 
call the Peierls - DeWitt bracket, has all the structural properties expected from a good 
Poisson bracket: bilinearity, antisymmetry, validity of the Jacobi identity and validity 
of the Leibniz rule with respect to plain and ordinary multiplication of functionals. 
This can be seen directly by noting that the first two properties and the Leibniz rule 
are obvious, while the Jacobi identity expresses the propagator identity for the causal 
Green function. But it is of course much simpler to argue that all these properties 
follow immediately from the above theorem, in combination with standard results of 
symplectic geometry. Moreover, the Peierls- DeWitt bracket trivially satisfies the fun- 
damental axiom of field theoretic locality: functionals localized in spacelike separated 
regions commute. All this suggests that the Peierls- DeWitt bracket is the correct clas- 
sical limit of the commutator of quantum field theory. Therefore, it ought to play an 
outstanding role in any attempt at quantizing classical field theories through algebraic 
methods, a popular example of which is deformation quantization. 

The basic complication inherent in the algebraic structure provided by the Peierls - 
DeWitt bracket is that it is inherently dynamical: the bracket between two functionals 
depends on the underlying dynamics. This could not be otherwise. In fact, it is the 
price to be paid for being able to extend the canonical commutation relations of classical 
field theory, representing a non-dynamical equal-time Poisson bracket, to a covariant 
Poisson bracket. The dynamical nature of covariant Poisson brackets is simplified (but 
still not trivial) for free field theories, where the equations of motion are linear, implying 
that the Jacobi operator 3[4>] an d its causal Green function G<j, do not depend on the 
background solution 0. 

Finally, we would like to remark that the main mathematical condition to be im- 
posed in order for the constructions presented here to work is that linearization of the 
field equations around any solution should provide a hyperbolic system of partial 
differential equations on M, for which existence and uniqueness of the Green functions 

and G^ can be guaranteed. There are various definitions of the concept of a hyper- 
bolic system that can be found in the literature, but the most appropriate one seems 
to be that of regular hyperbolicity, proposed by Christodoulou [28-30] in the context 
of Lagrangian systems, according to which the matrix 



u»u v 



8 2 L 
dql dqi 



should (in our sign convention for the metric tensor) be positive definite for timelike 
vectors u and negative definite for spacelike vectors u: a typical example is provided 
by the sigma model as discussed at the beginning of this subsection. What is missing 
is to translate this condition into the Hamiltonian formalism and to compare it with 
other definitions of hyperbolicity for first order systems, such as the traditional one of 
Friedrichs. 
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4 Conclusions and Outlook 



The approach to the formulation of geometric field theory adopted in this paper closely 
follows the spirit of Ref. [8], in the sense of emphasizing the importance of combining 
techniques from multisymplectic geometry with a functional approach. The main novel- 
ties are (a) the systematic extension from a Lagrangian to a Hamiltonian point of view, 
preparing the ground for the treatment of field theories which have a phase space 
but no configuration space (or better, a phase bundle but no configuration bundle), 
(b) a clearcut distinction between ordinary and extended multiphase space, which is 
necessary for a correct definition of the concept of the covariant Hamiltonian and (c) the 
use of the causal Green function for the linearized operator as the main tool for finding 
an explicit formula for the Hamiltonian vector field associated with a given functional on 
covariant phase space. This explicit formula, together with the resulting identification 
of the canonical Poisson bracket derived from the standard symplectic form on covariant 
phase space with the Peierls - DeWitt bracket of classical field theory, are the central 
results of this paper. 

An interesting question that arises naturally concerns the relation between the 
Peierls - DeWitt bracket as constructed here with other proposals for Poisson brackets 
in multisymplectic geometry. In general the latter just apply to certain special classes 
of functionals. One such class is obtained by using fields to pull differential forms / on 
extended multiphase space back to space-time and then integrate over submanifolds E 
of the corresponding dimension. Explicitly, in the Lagrangian framework, 

F[0] = / (EGofoty))*/, (115) 

whereas in the Hamiltonian framework, 

F[<t>] = [ (WofoTr))*/. (116) 

For the particular case of differential forms / of degree n — 1 and Cauchy hypersurfaces 
as integration domains S, this kind of functional was already considered in the 1970's 
under the name "local observable" [7] (though on ordinary instead of extended multi- 
phase space), but it was soon noticed that due to additional restrictions imposed on the 
forms / allowed in the construction, the class of functionals so defined is way too small 
to be of much use for purposes such as quantization. One of these restrictions is that / 
should be what is nowadays called a Hamiltonian form [32]. Briefly, an (n — l)-form / 
on J X ®E is said to be a Hamiltonian form if there exists a (necessarily unique) vector 
field Xj- on J l ®E, called the Hamiltonian vector field associated with /, such that 

i x u = df . (117) 

What seems to have motivated this restriction is the possibility to use the multi- 
symplectic analogue of the standard definition (|73|) of Poisson brackets in mechanics 
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for defining the Poisson bracket between the corresponding functional [8]. However, 
it turns out that, in contrast to mechanics where / is simply a function, the validity 
of eq. (jll7)) imposes strong constraints not only on the vector field Xj- but also on the 
form /; in particular, it restricts the coefficients both of Xj and of / in adapted local 
coordinates to be affine functions of the multimomentum variables p% and the energy 
variable p [33]. (See Refs [34,35] for a detailed analysis of the general situation encoun- 
tered when dealing with the same question for forms of arbitrary degree.) This implies 
that the class of functionals F derived from Hamiltonian (n — l)-forms / according to 
eqs (I115j) and/or (jllfjjl does not close under ordinary multiplication of functionals. 

Fortunately, using the Peierls- DeWitt bracket between functionals, we may dis- 
pense with the restriction to Hamiltonian forms. In fact, this line of reasoning was 
already followed by the authors of Ref. [8], where both the symplectic form on the 
solution space and the corresponding Poisson bracket between functionals on the so- 
lution space, with all its structurally desirable properties, are introduced explicitly. 
What remained unnoticed at the time was that this bracket is just the Peierls - DeWitt 
bracket of physics and that incorporating the theory of "local observables" into this 
general framework results in the transformation of a definition, as given in Ref. [7], into 
a theorem which, in modern language, states that the Peierls- DeWitt bracket {F, G} 
between two functionals F and G derived from Hamiltonian (n — l)-forms / and g, 
respectively, is the functional derived from the Hamiltonian (n — l)-form {/, g}. An 
explicit proof, based on the classification of Hamiltonian vector fields and Hamiltonian 
{n — l)-forms that follows from the results of Ref. [35], has been given recently [36]; 
details will be published elsewhere. 

Of course, there is a priori no reason for restricting this kind of investigation to forms 
of degree n — 1, since physics is full of functionals that are localized on submanifolds of 
space-time of other dimensions, such as: values of observable fields at space-time points 
(dimension 0), Wilson loops (traces of parallel transport operators around loops) in 
gauge theories (dimension 1), etc.. This problem is presently under investigation. 



Appendix: Affine Spaces and Duality 

In this appendix, we collect some basic facts of linear algebra for affine spaces which 
are needed in this paper but which do not seem to be readily available in the literature. 

A (nonempty) set A is said to be an affine space modelled on a vector space V if 
there is given a map 

+ : A , X \ - I (US) 
(a,v) i — ► a + v 

satisfying the following two conditions: 

• a + (u + v) = (a + u) + v for all a G A and all u, v e V . 
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• Given a, b G A, there exists a unique v G V such that a = 6 + t> . 



Elements of A are called points and elements of V are called vectors, so the map (jll8j) 
can be viewed as a transitive and fixed point free action of V (as an Abelian group) on A, 
associating to any point and any vector a new point called their sum. Correspondingly, 
the vector v whose uniqueness and existence is postulated in the second condition is 
often denoted by a — b and called the difference of the points a and b. 

For every affine space A, the vector space on which it is modelled is determined 
uniquely up to isomorphism and will usually be denoted by A. 

A map / : A — > B between affine spaces A and B is said to be affine if there exists 
a point a G A such that the map f a :A^B defined by 

f a (v) = f(a + v)-f(a) (119) 

—* — * — * 

is linear, that is, f a G L(A, B). It is easily seen that this condition does not depend on 
the choice of the reference point: in fact, if the map f a is linear for some choice of a, 
then the maps f a > are all equal as a' varies through A, so it makes sense to speak of 
the linear part f of an affine map /. Denoting the set of all affine maps from A to B 
by A(A, B), we thus have a projection 

I : A(A,B) — > L(AB) ^ (12Q) 

This construction is particularly important in the special case where B is itself a vector 
space, rather than just an affine space. Given an affine space A and a vector space W, 
the set A(A, W) of affine maps from A to W is easily seen to be a vector space: in fact 
it is simply a linear subspace of the vector space Map(A, W) of all maps from A to W. 
Moreover, the projection 

/ : A(A, W) — L(A,W) . . 

f — ► / 1 ] 

is a linear map whose kernel consists of the constant maps from A to W. Identifying 
these with the elements of W itself, we obtain a natural isomorphism 

A(A,W)/W = L(A,W) , (122) 

or equivalently, an exact sequence of vector spaces, as follows: 

— > W — ► A(A,W) L{A,W) — ► 0. (123) 

In the general case, one shows that given two affine spaces A and B, the set A(A, B) of 

► 

affine maps from A to B is again an affine space, such that A(A,B) = A(A,B), and 
that the projection ()12(JJ) is an affine map. 
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Concerning dimensions, we may choose a reference point o in A which provides not 
only an isomorphism between A and A but also a splitting of the exact sequence ()123j) 
and hence an isomorphism between A(A, W) and W © L(A, W), to show that 

dim A(A,W) = dimW + dim L(A,W) . (124) 

Choosing W to be the real line R, we obtain the affine dual A* of an affine space A: 

A* = A(A,R) . (125) 

Observe that this is not only an affine space but even a vector space which, according 
to eq. (|123|) . is a one- dimensional extension of the linear dual A* of the model space A 
by M, that is, we have the following exact sequence of vector spaces: 

— > K — > A* -U A* — > . (126) 

In particular, according to eq. ()124j) . its dimension equals 1 plus the dimension of the 
original affine space: 

dim A* = dimA + 1 . (127) 

More generally, we may replace the real line R by a (fixed but arbitrary) one-dimensional 
real vector space R (which is of course isomorphic but in general not canonically iso- 
morphic to M.) to define the twisted affine dual A® of an affine space A: 

A® = A(A, R) . (128) 

Again, this is not only an affine space but even a vector space which, according to 
eq. (J123)) . is a one-dimensional extension of the linear dual A* of the model space A 
by R, that is, we have the following exact sequence of vector spaces: 

— > R — > A® -U A® — ► . (129) 

Obviously, the dimension is unchanged: 

dim A® = dimA + 1 . (130) 

Moreover, we have the following canonical isomorphism of vector spaces 

A® ^ A" <g> R , (131) 

and more generally, for any vector space W, 

A(A,W) = A* ®W . (132) 

Regarding the splittings of the exact sequence (jl29|) . we note the following fact which is 
used in the construction of the inverse Legendre transformation: these splittings form 
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— * — * 

an affine space modelled on the bidual A** of A, which in finite dimensions can be 

— * 

identified with A itself. 

The concept of duality applies not only to spaces but also to maps between spaces: 
given an affine map / : A — > B between affine spaces A and B, the formula 

(/*(6*))(a) = 6*(/(a)) for b*eB\aeA (133) 

yields a linear map /* : B* — > A* between their affine duals 5* and A*. As a result, the 
operation of taking the affine dual can be regarded as a (contravariant) functor from the 
category of affine spaces to the category of vector spaces. This functor is compatible 
with the usual (contravariant) functor of taking linear duals within the category of 
vector spaces in the sense that the following diagram commutes: 

B* -A A* 

i I (134) 

B* -L-> A* 



Finally, we also need the construction of quotients in the affine category. These are 
defined by dividing out not affine subspaces but rather linear subspaces of the model 
space. In fact, given an affine space A and a linear subspace V of its model space A, we 
can declare two points a and a' of A to be equivalent modulo V if a — a' E V . Obviously, 
this relation is reflexive, symmetric and transitive, and hence is an equivalence relation 
dividing A into equivalence classes; the set of equivalence classes is as usual denoted by 
A/V. It is then easy to see that there is a unique affine structure on A/V turning A/V 

into an affine space such that A/V = A/V and such that the natural projection 

9 ■ A ^ A /Y (135) 
a i — > \a\ 

is an affine map. Moreover, this construction satisfies the standard factorization prop- 

— * — * 

erty: given two affine spaces A and B, two linear subspaces V of A and W of B and an 
affine map / : A — > B whose linear part / : A — ► B maps V into W, there exists a 
unique affine map [/] : A/V — > B/W such that the diagram 

/ 

A — > B 



P 



P 



(136) 



A/V — ► B/W 



commutes. 
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Concluding this appendix, we would like to point out that all the concepts intro- 
duced above can be extended naturally from the purely algebraic setting to that of 
fiber bundles. For example, affine bundles are fiber bundles modelled on an affine space 
whose transition functions (with respect to a suitably chosen atlas) are affine maps. 
Moreover, functors such as the affine dual or the construction of quotients are smooth 
(see [31] for a definition of the concept of smooth functors in a similar context) and 
therefore extend naturally to bundles (over a fixed base manifold M). This means that 
any affine bundle A over M has a naturally affine dual, which is a vector bundle A* 
over M, and that given any vector subbundle V of the difference vector bundle A of an 
affine bundle A over M, we can form the quotient affine bundle A/V over M. 
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